Reconfigurable Gradient Index using VO2 Memory Metamaterials 
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We demonstrate tuning of a metamaterial device that incorporates a form of spatial gradient control. Electrical 
tuning of the metamaterial is achieved through a vanadium dioxide layer which interacts with an array of 
split ring resonators. We achieved a spatial gradient in the magnitude of permittivity, writeable using a single 
transient electrical pulse. This induced gradient in our device is observed on spatial scales on the order of 
one wavelength at 1 THz. Thus we show the viability of elements for use in future devices with potential 
applications in beamforming and communications. 



Metamaterials have progressed from academic 
curiosities 1 2 to candidates for real- world applications P 4 ^ 
Emerging metamaterial applications range from radio 
frequency (RF)P communications to millimeter radar. 6 
One key technique which promises to further the ap- 
plicability of metamaterials is tunability. Operational 
frequency tuning has been presented as one solution 
to the narrow bandwidth often present in metamate- 
rial devices Frequency- agile metamaterial designs 
have been demonstrated across a wide spectrum from 
microwave^ to near- visible frequencies.^ 

To date, tuning has generally been homogeneously 
implemented across the entire device as a whole. De- 
veloping techniques for controllable spatially variable 
tuning will present the possibility of devices with a 
reconfigurable Gradient Index of Refraction (GRIN). 
GRIN devices have already proven an attractive area for 
metamaterialsp^as the metamaterial design process nat- 
urally allows for the control needed to fabricate GRIN 
structures. Additionally, use of spatially nonuniform tun- 
ing can leverage the narrow bandwidth of metamaterials 
in a unique way. For narrow-band operation, minor ad- 
justments in the resonance frequency of a metamaterial 
can translate to large changes in the index of refraction at 
that frequency. Overall, spatial control of resonance tun- 
ing allows for post-fabrication modification of the index 
of refraction and therefore the creation of a reconfigurable 
gradient. 

In this work, we demonstrate a spatially reconfig- 
urable THz hybrid metamaterial with vanadium dioxide 
(VO2) and split ring resonators (SRRs) as constituent 
elements. The SRR has been the "fruit fly" of metama- 
terials research, allowing for convenient implementation 
of optical characteristics which are unattainable without 
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metamaterials P3I2113I Our device is composed of an array 
of 100 nm thick gold SRRs (dimensions shown in fig- 
ure [5| lithographically fabricated on 90 nm thick VO2 
grown on a sapphire substrate. VO2 undergoes an insu- 
lator to metal transition^ which can be triggered ther- 
mally electrical!}^ or optically.^ The phase transition is 
hysteretic, and therefore, changes in the conductivity of 
VO2 generally persist, provided the device temperature 
is maintained. Hybrid metamaterial- VO2 devices benefit 
from this memory,^ and from the large tuning dynamic 
range achievable with VO2. Persistent tuning uses this 
memory to eliminate the need for continuous stimulation 
or repeated excitation of the device. Our measurements 
demonstrate the ability to perform persistent tuning of 
a hybrid SRR-VO2 device with a spatial configuration, 
writing a gradient in the permittivity of the hybrid meta- 
material. 

To interrogate our metamaterial, we performed in- 
frared transmission spectroscopy using a home-built 
broadband IR microscope. Linearly polarized light from 
a mercury lamp was focused via reflecting optics to a 
normally-incident spot ~100/im in diameter. In far- 
infrared, this spot facilitates nearly diffraction-limited 
measurements of the metamaterial local response. The 
sample was mounted on a translator allowing for linear 
movement across the entire device. Integration of a tem- 
perature stage with proportional-integral-derivative con- 
troller (PID) feedback allowed us to control device tem- 
perature to better than ±0.05 K. For these experiments, 
we kept the device at a constant temperature of 339.5 K 
(just below the temperature of the insulator-metal tran- 
sition in VO2) to allow for maximal transition during the 
voltage application discussed below. 

Prior to the application of any current pulse, cali- 
bration transmission spectra were collected at positions 
across the sample at 339.5 K. In these spectra, the peak 
position of 1- T(u) (transmission) indicates the resonance 
frequency of the hybrid metamaterial. In our setup, the 
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device was oriented to electrically excite the lowest or- 
der RLC-oscillation mode of the SRR. The increasing 
permittivity of VO2 as its phase transition progresses in- 
creases the capacitance of each SRR. This in turn reduces 
the effective time constant of the RLC-oscillation mode 
thereby decreasing the resonance frequency of the hybrid 
SRR-VO2 device. Additionally, losses in the metallic re- 
gions of VO2 cause damping of the resonance. These 
calibration spectra are nearly identical with uniform reso- 
nance positions, indicating the homogeneous (on a meso- 
scopic effective medium scale) initial state of our device. 

To electrically write a spatial gradient to the device, 
we applied a square voltage pulse of 175 V x 2.5 s across 
the contact points as shown in Figure [I] (top inset). This 
voltage pulse did not significantly affect average tempera- 
ture of the device as the total input energy is small com- 
pared to the thermal inertia of the sapphire substrate. 
However, as current flowed almost exclusively through 
the VO2-SRR layer, local transient heating of this layer 
occurred. While the dominant switching mechanism in 
this device is thermal, there may be an additional contri- 
bution due to electrical switching from the applied volt- 
age. Inquiry into the microscopic me chani sms of switch- 
ing is an active area of VO2 researcrP^DH DU t i s beyond 
the scope of our work, which focuses on the utility and 
applications of switching effects. As VO2 undergoes its 
transition from insulator to metalf^we expect the peak 
of the resonance to redshift and become damped. 19 Such 
transient heating leave s a persistent change due to the 
hysteresis in VC^P^^ 

The placement of electrical contacts was designed to in- 
duce a gradient in the current density passing through the 
VO2-SRR layer. Using our translation stage, we probed 
8 locations on the sample, ranging from directly between 
the contacts (0 mm) to +3.5 mm above the contact path 
as shown in Figure [l] (top inset). In Figure T] we plot 
transmission spectra at these 8 positions. The redshift 
of the absorption resonance is more prominent for posi- 
tions in close proximity to the contacts. This is consis- 
tent with the idea that the (percolative) phase transition 
in VO2 has progressed further near the contacts. The 
maximum resonance redshift is about 10% of the center 
frequency, and nearly 50% of that obtained from the en- 
tire insulator-metal transition!^ Through a comparison 
of the measured resonance frequencies to those in a previ- 
ous work using the same device,^ we estimated an effec- 
tive temperature gradient, during and shortly after the 
current pulse, ranging from ~342 K at mm, to ~340 K 
at 3.5 mm. Additionally, this experiment demonstrates 
that metamaterials can be probed and meaningfully mod- 
ified on spatial scales on the order of the wavelength. The 
total spatial extent of our observed gradient is 3.5 mm, 
with 8 data points over this distance. Thus, these dis- 
tinct spectra are recorded with a spacing of just under 
1.5A (~ 300/im at ITHz). 

To better quantify the potential of such gradients, we 
retrieved the optical parameters from the transmission 
curves by fitting. Permittivity gives a better visualization 
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FIG. 1. (Color online) Transmission measurements taken 
along the hybrid metamaterial at 339.5 K. mm indicates 
a position midway between the contacts. (Bottom inset) 
Measured transmission curves along with fits. (Top inset) 
Schematic of device indicating contact geometry and mea- 
surement locations. 



of the gradient in material response than transmission, 
especially in a single-layer metamaterial where overall 
transmission is generally high. The retrieval procedure is 
a two-oscillator fitting using Lorentzian o scillat ors e and 
pi weighted by spatial-dispersion cofactorj 2Q | 21> to give a 
close approximation to the homogenized effective permit- 
tivity, e e ff, and permeability, /^ e ff of the metamaterial 
layer: 
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where is the phase advance across one cell, and Sb=l 
for electric resonators or St=-l for magnetic resonators. 
The resonance frequency, strength, and damping con- 
stant of two oscillators 22 underlying the effective optical 
parameters were modified to obtain good agreement be- 
tween measured and fitted spectra. While better agree- 
ment may have been obtained through the use of a large 
number of oscillators, this procedure increases numerical 
complexity, and can often obscure insight into the un- 
derlying physical processes involved in the material re- 
sponse, here the magnetic (u; ) and electric (u>i) modes 
of the SRRpSl These oscillator parameters were then used 
to determine permittivity and permeability of our de- 
vice. Resultant fits for the most shifted (0 mm) and least 
shifted (3.5 mm) transmission lines are shown in Figure 
[l] (bottom inset). The fits capture the most prominent 
features of the transmission spectra. They also model the 
pulse-induced redshift and damping of the resonance. 
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FIG. 2. (Color online) Extracted values for real (a) and imag- 
inary (b) permittivity. Position of maximum permittivity in- 
dicated by dashed lines. (Inset) Schematic of SRR unit cell. 



The real and imaginary permittivity extracted from 
these fits are shown in Figure |2j These figures illustrate 
the range of complex permittivity values simultaneously 
present at different locations in the device. By inter- 
polating between the 8 points where the spectra were 
obtained, we extracted a spatial map of the metamate- 
rial permittivity response, shown in Figure [3] (a). At 
locations nearer to the contacts, the resonance becomes 
increasingly redshifted due to the increasing number and 
density of metallic puddles in VO2. The overall frequency 
tuning range of our device compares favorabl y with many 
techniques for frequency- agile metamaterialsP^'^' At the 
initial resonance frequency (~1.23 THz), the overall per- 
cent change in real permittivity is 70% across the sample, 
enough to satisfy many potential GRIN device applica- 
tions. 

The shift in resonance position in our transmission 
spectra is a result of the application of a transient voltage 
pulse as described previously. Current flowing through 
the device heats the VO2 with power dissipated given as 
P = I 2 R. We expect that areas of the device nearer 
to the contacts were exposed to a larger current density, 
and as a result the VO2 will have progressed through a 
greater amount of its phase transition. To model this 
situation, we performed a simulation of current flow and 
power dissipation in our device using COMSOL, a finite 
element simulation program. The results of this simula- 
tion are shown in Figure [3] (b) . As expected, the great- 
est power dissipation occurs in areas near the contacts. 
These results fit with our picture of a spatial gradient of 
permittivity induced by inhomogeneous excitation of the 
VO2 phase transition. 

Our SRR-VO2 device is the simplest hybrid meta- 
material structure with which to explore advantages of 
spatial gradient tuning. With it, we have shown the 
ability to induce and spectroscopically probe a persistent 
spatial gradient within a previously uniform device. Use 
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FIG. 3. (Color online) a) Spatial map of real permittivity 
varying with vertical position and frequency. Resonance fre- 
quency locations indicated by colored points, b) Simulation 
of current flow through device. Color indicates power density 
(W/m 3 ) and current flow is shown in the vector field (A/m 2 ). 
Black rectangular outline indicates the location of the SRR 
array and VO2. Contact locations and shapes indicated by 
outlines at the bottom of the device. 



of finer, or even pixel- by-pixel, control over switching 
would allow for creation of more precise gradients. Such 
pixel-level switching nearly necessitates the existence of 
memory within the device as the difficulty of repeated 
or sustained excitation would make detailed switching 
unfeasible. Future flexibility in the optical parameters of 
metamaterials as a result of tuning occurring on the pixel 
level will enable the creation of dynamic hybrid metama- 
terial structures suitable for a wide range of applications. 

The authors wish to acknowledge funding for this 
project from AFOSR and the creative research project 
of ETRI. 
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